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Abstract. The Gribov ambiguity problem is studied for compact lattice 
QED within the Lorentz gauge. In the Coulomb phase, Gribov copies are 
mainly caused by double Dirac sheets and zero-momentum modes of the 
gauge fields. Removing them by (non-) periodic gauge transformations al- 
lows to reach the absolute extremum of the Lorentz gauge functional. For 
standard Lorentz gauge fixing the Wilson fermion correlator turns out to 
be strongly effected by the zero-momentum modes. A reliable fermion mass 
estimate requires the proper treatment of these modes. 



1. Introduction 

Most of the applications of lattice gauge theories are based on and are 
employing their manifest gauge invariance. However, in order to get a bet- 
ter understanding of the structure of the lattice theory itself and to inter- 
prete correctly results obtained in Monte Carlo simulations, it is instructive 
to compare also gauge variant quantities such as gauge and fermion field 
propagators with corresponding analytic perturbative results. In this re- 
spect, compact lattice QED within the Coulomb phase serves as a very 
useful 'test ground'. In the weak coupling limit this theory is supposed to 
describe non-interacting massless photons. 
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In order to fix the gauge the Lorentz (or Landau) gauge condition is 
normany apphed. For non-Abehan gauge theories there is no unique solu- 
tion, i.e. so-called Gribov copies occur [1]. Within continuum QED such a 
problem arises, too, if the theory is defined on a torus [2]. The lattice dis- 
cretization may cause additional problems. Indeed, various lattice studies 
[3 - 11] have revealed nontrivial effects. The standard Lorentz (or Landau) 
gauge fixing procedure leads to a r-dependence of the non-zero-momcntum 
transverse photon correlator inconsistent with the expected zero-mass be- 
havior [3]. Numerical [4, 5, 9] and analytical [7] studies have shown that 
there is a connection between 'bad' gauge (or Gribov) copies and the ap- 
pearance of periodically closed double Dirac sheets (DDS). The removal 
of DDS by appropriate gauge transformations restores the correct pertur- 
bative behavior of the photon correlator at non-zero momentum, but it 
does not completely resolve the Gribov ambiguity problem. Gribov copies 
related to other local extrema of the gauge functional and connected with 
zero-momentum modes (ZMM) of the gauge fields still appear. They 'dam- 
age' gauge dependent obscrvables such as the zero-momentum gauge field 
correlator [7, 9, 8] and the fermion propagator [10, 11], respectively. 

There is a special Lorentz gauge, for which both the double Dirac sheets 
and the zero-momentum modes can be removed from the gauge fields . We 
call it zero-momentum Lorentz gauge (ZML) [8]. It allows to reach the 
global extremum of the Lorentz gauge functional in almost 100% of the 
cases. In comparison with the standard Lorentz gauge procedure (LG) it 
demonstrates very clearly the strong effects caused by the zero-momentum 
modes. 

In the given talk we are going to review the results of [8, 11] with 
special emphasis on the question, how Gribov copies influence the Wilson- 
fermion propagator within the Coulomb phase of quenched QED. We want 
to show that a reliable estimate of the (renormalizcd) fermion mass requires 
either the removal of the zero-momentum modes or their proper perturba- 
tive treatment, when comparing the numerical results with analytic expres- 
sions. By employing the ZML-gauge we shall estimate the fermion mass in 
agreement with standard perturbation theory. 

2. The Action and the Observables 

We consider 4d compact QED in the quenched approximation on a finite 
lattice (y = Ng x Nt). The pure gauge part of the standard Wilson action 
[12] reads 



Sg = P (1 ~ cos6'j,,^,,) , 

x,ij,<iy 



(1) 
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with the plaquette angle 9x,,xy = Ox,ii + Ox+fi,y - 6x+i',^i - Ox,u related to the 
link variables 9x,ix G (— 7r,7r]. /? = l/eg is the inverse bare coupling. The 
lattice spacing is put a = 1, for simplicity. 
The fermion part is given by 

SF = Y.'^x^xy{e)^y, M = 1 - «D, (2) 

x,y 

4 

^xy = 5Z { ^y,x+p, + ^r-/i,//-^/i"^2/,3;-/i} j 

where = 1 ± 7^ and Ux,i^ = e'^^ ''. The hopping-parameter k is related 
to the bare mass uiq by At = 1/(8 + 2mo). 

In quenched QED the observables have to be averaged with respect 
to the gauge field {9} with the weight exp(— 5g) . We imply periodic 
boundary conditions (b.c.) except for the fermion fields, which we choose 
to be anti-periodic in the X4 = r-direction. 

The first gauge variant observable we are going to discuss is the trans- 
verse photon correlator at non-zero momentum 

T^\p;t) = ($(p;r)$*(p;0)) , (3) 

X 

with { jj, = 1,3, p = {0,p, 0) ). The second one is the fermion propagator. 
For a given gauge field {9} we have 

x,X4 y 

In the following we shall restrict ourselves to the vectorial part 

ry(T) = ^ReIV(74r(T)), (5) 

with the trace taken with respect to the spinor indices. For the b.c. men- 
tioned above, (ry(T)) is an even function of r — Nt/2. 

Lateron, we shall compare the expectation value ( Ty ) with the result 
of a simple approximation, which takes only constant gauge field modes into 
account. The correlator in a uniform background 9x,fj, = 4'ixi ~ tt < 0^ < tt, 
/X = 1, • • • , 4 can be represented as 

[gr ^ g2jVt-r] cos(04r) + + g^''^] COs[4>i{Nt - t)] 

^ 1 + f2iVi + 2£^^ COs((/)4iVt) ' 
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where 



£ = l + 



M^ + K? _ ^/M'^ + K?^J{M + 2)^^ + K:^ 

2(1 2{l + M) 



3 



3 



= mo + ^ (1 - COS (t)i) , 
1=1 



^sin^ 4>i, mo> 0. 
1=1 



For (pfj, = 0, fi = 1, ■ ■ ■ ,4 the free fermion correlator for finite lattice size 
[13] is reproduced. 

3. Lorentz Gauge Fixing 

In numerical simulations the Lorentz gauge is fixed by iteratively maximiz- 
ing the gauge functional 



The algorithm is called standard Lorentz gauge fixing (LG), if it consists 
only of local maximization and overrelaxation steps [14] with respect to 
gauge transformations periodic in space-time. The standard procedure gets 
normally stuck into local maxima of the gauge functional (7) (gauge copies). 
It has been argued that the Gribov problem has to be solved by search- 
ing for the global maximum providing the best gauge copy [15]. In [8] we 
have shown that in order to reach the global maximum we have necessarily 
to suppress both the double Dirac sheets (DDS) and the zero-momentum 
modes (ZMM) in the gauge fields. Let us explain this in more detail. 

DDS can be identified as follows. The plaquette angle 9x,^ii' is decom- 
posed into the gauge invariant (electro-) magnetic flux 9x,^j,i/ G (~^5 ^] a^nd 
the discrete gauge-dependent contribution 2iTnx^fj,u, n^^f^u = 0,±1,±2 [16]. 
The latter represents a Dirac string passing through the given plaquette if 
nx,^^ = ±1 {Dirac plaquette). A set of Dirac plaquettes providing a world 
sheet of a Dirac string on the dual lattice is called Dirac sheet. DDS consist 
of two sheets with opposite flux orientation extending over the whole lattice 
and closing themselves by the periodic b.c. They can easily be identified by 
counting the total number of Dirac plaquettes Nj^p^ for each choice (/x; u). 
The necessary condition for the occurence of DDS is that at least for one 
of the six possibilities (/x; v) holds 




with respect to the (local) gauge transformations 

Ux^, — ^ AxUx^,Al^. ; = exp{iJ^^} e U{1) . 



(8) 



(9) 
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Figure 1. Time history of Fmax (a) and Nop (b) at /3 = 1.1 on the 12 • 6^ lattice in 
the standard Lorentz gauge. 

DDS can be removed by periodic gauge transformations. 
The ZMM of the gauge field 

= 7 E^-,^ (10) 

X 

do not contribute to the pure gauge field action either. For gauge configu- 
rations representing small fluctuations around constant modes it is easy to 
see, that the global maximum of the functional (7) requires (f)^,^ = . The 
latter condition can be achieved by non-periodic gauge transformations 

Ox,u. = Cn + Gx,n mod 27r, G (-tt, tt]. (11) 

We realize a proper gauge fixing procedure as proposed in [8]. Successive 
Lorentz gauge iteration steps are always followed by non-periodic gauge 
transformations suppressing the ZMM. Additionally we check, whether the 
gauge fields contain yet DDS. The latter can be excluded by repeating 
the procedure with initial random gauges. We call the combined procedure 
zero-momentum Lorentz gauge (ZML gauge) . It yields the global maximum 
of the gauge functional with very high accuracy. 
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Figure 2. Time history of Fmax (a) and Nnp (b) at = 1.1 on the 12 • 6^ lattice in 
ZML gauge. 



In Figures 1 and 2 we show, how the achieved values of the gauge func- 
tional (7) are correlated with the occurence of DDS visible as sharp peaks 
in the number of Dirac plaquettes. Whereas for LG strong fluctuations oc- 
cur, they disappear after ZML gauge. The few DDS seen in Fig. 2 are easily 
removed by restarting the procedure with random initial gauges. Random 
gauges can also be used in order to convince oneself that the ZML gauge 
prescription leads to the global maximum of the gauge functional in more 
than 99% of the cases. 

4. Results 

First let us convince ourselves that the removal of the above mentioned 
gauge copies leads to the correct behaviour of the transverse photon propa- 
gator. In Fig. 3 we show the normalized correlator T^\p; T)/T^{p; 0) for 
lowest non-vanishing momentum and for different Lorentz gauge prescrip- 
tions. For the standard one (LG) we see a clear deviation from the expected 
perturbative zero- mass result. We show also the result obtained with an 
axial Lorentz gauge (ALG) using an initial maximal-tree axial gauge con- 
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Figure 3. Transverse propagator at (3 = 1.1 on the 12 ■ 6^ lattice in three different 
gauges. 

dition [17], which provides a 'unique' prescription. The latter gauge fixing 
prescription turns out to be even worse! On the other hand the ZML gauge 
provides an excellent agreement with the perturbative result. In fact, as we 
convinced ourselves earlier, it is already sufficient to remove the DDS gauge 
copies in order to reach this agreement [4]. The given observations do not 
change, when /3 and/or the lattice size are increased considerably [9]. 

In the following we want to concentrate on the pure effect of the ZMM. 
Therefore, we compare the ZML gauge with a version of the standard 
Lorentz gauge, where the DDS are removed and the ZMM arc left. We 
shall abbreviate the latter version also by LG. For both the LG and ZML 
gauges we have computed the averaged fermion correlator employing the 
conjugate gradient method and point-like sources. In the upper part of Fig. 
4 we have plotted Tv{t) (normalized to unity at r = 1). The situation 
seen is typical for a wide range of parameter values within the Coulomb 
phase. Obviously, there is a strong dependence of the fermion propagator 
on the gauge fixing procedure resulting in the presence or absence of ZMM. 
The masses to be extracted seem to have different values. Let us determine 
the effective mass me//(r) in accordance with 

(r(T + 1; 9))e ^ cosh[E{T){Nt/2 - r - 1)] 
(r(r;0)), cosh[E{T){Nt/2-T)] ^ ^ 

where E{t) = ln{meff{r) + 1). See the lower part of Fig. 4. In the LG 
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Figure 4- The fermion vector propagator (a) and the effective mass (b) at /3 = 2 and 
K = .122 on a 12 X 6^ lattice for LG and ZML gauges as explained in the text. 



case no plateau is visible, whereas the ZML case provides a very stable 
one. Thus, only the ZML gauge yields a reliable mass estimate, whereas 
the standard method to fix the Lorentz gauge obviously fails. 

To get deeper insight into the effect of ZMM for the LG case (with DDS 
suppressed) we measure the probability distributions P{4>) for the space- 
and timc-likc components of ZMM according to Eq. (10). The distributions 
turn out to be flat up to an effective cutoff at |(/)^| ~ and to be 

widely independent of (3. In accordance with Eq. (6) we compute the fermion 
propagator for constant modes in the LG case and average 

{Tvir; 0))<^ = I m Pi4>) Ty (r; </>). (13) 

The results for several parameter sets are presented in Fig. 5 together with 
the corresponding free (i.e. zero-background) propagator. We see clearly 
that the constant mode contributions strongly change the behavior of the 
fermion propagator and, naively speaking, produce a larger mass. 
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Figure 5. Free fermion propagator (dashed line) and averaged constant-mode propagator 
in the LG case (full line) for (3 = 2,10, k = .120, .124, lattice size 12 x 6^, 16 x 8^. 

Finally, in Fig. 6 we present the fermion mass extracted from the vector 
fermion propagator within the ZML gauge for /? = 2.0 and various k- values. 
We see a nice linear behaviour from which by extrapolating to zero mass 
(solid line) we estimate the critical value Kc = 0.1307 ± 0.0001. 

5. Conclusions 

We have studied the effect of different gauge copies of the gauge field on 
gauge dependent correlators, in particular on the Wilson fermion propaga- 
tor. We have convinced ourselves that the standard Lorentz gauge fixing 
prescription to maximize the functional (7) provides gauge copies with DDS 
and ZMM. These modes disturb the photon and the fermion correlator in 
comparison with perturbation theory and consequently spoil the (effective) 
mass estimate. A Lorentz gauge employing non-periodic gauge transfor- 
mations in order to suppress the ZMM - additionally to DDS - (the ZML 
gauge) allows to reach the global maximum of the Lorentz gauge functional. 
Furthermore, it provides a reliable fermion mass determination, at least, if 
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Figure 6. Fermion mass as a function of inverse n obtained within the ZML gauge 
for (5 — 2.0onal2x6^ lattice. The solid line represents a linear fit providing 
Kc{/3) = 0.1307 ±0.0001 . 



K is chosen not too close to the chiral critical hne Kc{P)- A computation of 
the fermion propagator with constant background gauge fields taken from 
the ZMM of the quantum fields demonstrates the disturbing effect of these 
modes very clearly. Moreover, it shows the effect to be independent of the 
bare coupling and not to disappear for large volumes. 

So far, we have studied the quenched approximation of U(l) lattice 
gauge theory. The gauge action (1) is invariant under non-periodic gauge 
transformations (11). Thus, we are allowed to use the ZML gauge for evalu- 
ating gauge dependent objects. Contrary to the gauge action, the fermionic 
part (2) does depend on the ZMM because of the (anti-) periodic boundary 
conditions. In this case another way of dealing with the Gribov problem 
has to be searched for. 

The problems we have discussed here for compact QED show that gauge 
fixing has to be carried out and to be interpreted with care. This lesson has 
to be taken into account also in lattice QCD when extracting masses from 
gauge variant gauge and fermion correlators, respectively. 
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